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Abstract 

After analyzing Dirac's equation, one can suggest that a well-known quantum-mechanical mo- 
mentum operator is associated with relativistic momentum, rather than with non-relativistic 
one. Consideration of relativistic energy and momentum expressions allows us to define the non- 
relativistic, relativistic and pseudo-relativistic (present in Schrodinger equation) kinetic energy 
operators. Consequences of kinetic energy operator's correction for spectra of basic quantum 
Hamiltonians are investigated. In some cases this correction can produce remarkable spectra mod- 
ifications. 
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I. INTRODUCTION 



The many-particle Schrodinger equation is the best tool for quantum systems low-energy 
phenomena investigation. The existing powerful approximate methods enable successful de- 
scription of such systems as electrons in solid state, molecules, atoms and even the atomic 
nuclei. The Hamiltonian is defined as a non-relativistic operator, hence relativistic modifica- 
tions are among the most important. However, the solution of this problem is very problem- 
atic due to the complex relativistic kinetic energy operator and the undefined many-particle 
Dirac equation. This problem has been solved only in weakly relativistic approximation for 
systems with Coulomb potential, such as an atom. The main points of this solution are 
based on Hamiltonian's corrections, saving the usual form of kinetic energy operator and 
introducing new potentials like spin-orbit, following from corresponding Dirac equation. For 
systems with strong interaction these Hamiltonian corrections are not sufficient in order to 
obtain acceptable result. New ways need to be investigated in order to solve this problem. 

Let us take a new look at relativistic modifications for basic quantum-mechanical Hamil- 
tonians describing particle, moving in spherically-symmetric external fields, created by deep 
well, harmonic oscillator and Coulomb potentials. These Hamiltonians consist of one-particle 
kinetic energy operator and external field potential operator, and can be presented as: 

where mc^ is particle's rest energy, given in GeV, he = 0.197 GeV fm. r is absolute value of 
it's radius- vector and r marks the set of intrinsic variables, such as mass, charge, spin, isospin 
and possible others. Such operators are the essential parts of many-particle Hamiltonian, 
hence relativistic corrections are necessary for further investigation of this problem. 



II. DIRAC EQUATION TRANSFORMATIONS 

Dirac equation is the basic for Hamiltonian's relativistic corrections: 

k=l / 

where ctk and (3 are Dirac matrices and wave-function is defined as set of four functions 
- components. Modified stationary equation for free particle with spin, equal h/2 can be 



written in Pauli matrices as system of two equations: 

(e - mc^) (f + ihc{(T-V)x = 0, 
{e + mc^) X + ifi'C (cr • V) (^9 = 0. 

Here e is relativistic particle's energy operator, and cr - spin operator, ip and x ^"^^ two- 
component spinors composing entire wave-function 



^ = r • (4) 




Let us define \I/ as eigenfunction of operator e corresponding to eigenvalue, equal e. Having 
in mind that commutator of operators e and V equals zero, the eigenfunction of e is also 
eigenfunction of momentum operator p = —iK^ with eigenvalue, equal p, one obtains that 



e 
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{m(?Y + {pcf = ['ymc'^Y . (5) 



This is the well-known expression for energy of the relativistic particle with momentum 
p = jniv, where Lorentz factor 

7={l-iv/cf)-'^\ (6) 

This leads to the conclusion that quantum-mechanical momentum operator p = —iK^ is 
associated with relativistic momentum p = 'ymv rather than with non-relativistic one Pq = 
mv. 

The non-relativistic momentum operator can be defined in a following way. From Eq. 
([5]), written in operators form, it follows that 

p2 = (7^ - 1) {mcf , (7) 

where the operator 

y = (l-(po/mcf)"\ (8) 

Defining this classical Lorentz factor as quantum-mechanical operator, gives us the pos- 
sibility to introduce non-relativistic momentum operator, non-relativistic, relativistic and 
pseudo-relativistic (present in Schrodinger equation) kinetic energy operators. As well as 
analyze relations between them and investigate spectrum modifications when different ki- 
netic energy operators are present in Hamiltonian. 



Applying the definition of tlie non-relativistic kinetic energy operator 



tn = Po/2m (9) 

one obtains tlie expression: 

P' = PoT' = 2mt„ (1 - 2tjmc^) ~' . (10) 

Let us introduce tlie operator present in Schrodinger equation and call it pseudo-relativistic 
kinetic energy operator: 

to-P>--igA. (U) 

It can be present in terms of t„ : 

to = t„(l-2t„/mc')"'. (12) 
From last expression it follows that 

2 

tn = 1^ (1 + {p/mcYy' = to (1 + 2to/mcy' , (13) 

or 



tn — to 



1 - 2to/mc2 + (2to/mc^) - ... . (14) 



Therefore, the non-relativistic kinetic energy operator tn is far more complex than pseudo- 
relativistic one to. 

Having kinetic energy definition in classical relativity t = e—mc^ in mind, one can present 
the system of Dirac equations ([3]) in the form of quantum-mechanical operators 

^ ^ (15) 
(t + 2mc^) X = c{(T ■ p) ip 

and after usual manipulations 

X= {t + 2mcY'c{(T-p)ip, (16) 

t(p = c{(T-p){t + 2mc^y^c{(7-p)ip, (17) 
can obtain the relativistic kinetic energy operator: 



'mc^)^ + (pc)^ — mc? . (18) 
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In terms of to it can be presented as 



+ 2to/mc2 - ij . (19) 



This can be written as a series expansion: 



t = tr 



1 - to/2mc2 + {to/mc^Y /2 - ... , (20) 



or 

'2k 



t-t V (-^)' i 



(to/mc^)'. (21) 



k 

The first term of all operators equals to, see f|T^ . fl20|) . hence at non-relativistic energies, 
when to ^ mc^, all operators have the same form. 

The expression of relativistic kinetic energy operator in terms of non-relativistic one is 

t/mc^ = (1 - 2t„/mc2)"^^^ - 1. (22) 

The commutators of all pairs of presented kinetic energy operators equal zero, all of them 
are differential operators of the same argument - radius vector of particle - hence relations 
between eigenvalues can be defined as relations between operators, and can be considered 
in terms of classic relativity. The eigenvalues, corresponding to the same eigenfunction, are 
distributed in a following way: 

tn<t< to. (23) 

Thus, the pseudo-relativistic operator produces the largest eigenvalue in comparison with 
other two for the same eigenfunction. All operators are positively defined, so it follows that at 
non-relativistic kinetic energy, corresponding to the maximum possible velocity v = c, when 
tn/mc^ = 1/2, both (relativistic ( l22l) and pseudo-relativistic ( TT2l) ) kinetic energy operator's 
eigenvalues approach, as necessary, the infinite values. 

Therefore, the Schrodinger equation is not completely non-relativistic. Some relativistic 
corrections have already been taken into account applying relativistic momentum operator 
instead of the non-relativistic one. In first approximation, i.e. at negligible in comparison 
with mc^ kinetic energy, all operators coincide. 

Dirac equation ([3]) can be transformed into Schrodinger equation form also for particle 
moving in external field, which scalar potential is v. It is necessary to introduce the ex- 
pression t = h — V, where h is the Hamiltonian instead of kinetic energy operator t. The 



equation (fTSi) can be rewritten in following form 

frr • \ I ld \ 

(24) 

Now let's assume that present wave-function is an eigenfunction of Hamiltonian 

(25) 



c (o- • p) 1 


II 




p) V— 2mc^ j 















l-l 











where e is corresponding eigenvalue. The system of equations now is: 

(e-v) = c (<T ■ p) X, ^26) 
e — V + 2mc^) x = c (cr • p) (p. 

As usual, the following operation eliminates the wave-function component x- The equation 
for component ip appears in form: 

ey? = WLf + c (cr ■ p) (e — v + 2mc^) c (cr ■ p) yj. (27) 

The commutator of potential energy operator v and momentum operator p does not equal 
zero. This makes the expression complex enough, so it can be successfully developed only 
at fulfilled condition e — v ^ 2mc^, known as weakly relativistic approach . After the 

simple manipulations it can be presented as 

h = ho+vr + vs + Vf), (28) 

where apart the standard Hamiltonian hg = to+v, the additional potentials appear. 

(e— V (r))^ , , 

zmc^ 

equals relativistic kinetic energy second order term f l20|) , converted to potential operator 
form. 

1 ^ he \ ^ 1) '^^ (30) 

2 \mc^ J r dr 

is the spin-orbit potential (here s and 1 are dimensionless spin and orbital momenta opera- 
tors) and 

^r> = U-^X^^ir) (31) 
8 \mc'^ J 

is the Darwin potential. 



This investigation shows that potential's weakly relativistic modification has to be con- 
sidered only as some approximation for complete kinetic energy conversion to relativistic 
form. 

Dirac equation for particle, moving in external field, transformation into Schrodinger 
equation form was very attractive years ago. Namely because kinetic energy expression 
is like one, present in Schrodinger equation, e is characteristic for Schrodinger equation 
energy of particle's bound state in given potential (equal binding energy with negative sign) 
and additional terms of a potential have a well-defined structure, appearing because of the 
nonzero spin and associated magnetic momentum. 

However nucleons or constituent quarks besides of spin have additional intrinsic degrees of 
freedom, such as isospin and color and the realistic potentials contain all operators, allowed 
by symmetry considerations, these weakly relativistic modifications of strong interaction 
potentials are not necessary. 

III. BASIC QUANTUM HAMILTONIAN'S RELATIVITY 

Let us now investigate bound-states spectrum of basic quantum-mechanical Hamiltonian 
(II]) when relativistic kinetic energy operator is applied instead of pseudo-relativistic one. 
The Schrodinger equation for these Hamiltonians is 



The spherically-symmetric potential provides the possibility to reduce this equation to or- 
dinary differential equation for radial function, defined in a following way: 



where (i^y?) is spherical harmonic and n = 1,2,... equals the number of bound state, 
counted from the lowest one. The radial equation after some transformations can be pre- 
sented as 



Here the function Uni (r) = rRni (r) is more acceptable than the radial function. Radial 
Schrodinger equation for this function has the form 




(32) 




(33) 




(34) 




(35) 
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The attracting points of this radial function are as follows: it is normalized by condition 



\uni {r)\'^dr = 1, 



(36) 



so \uni {r)f equals the radial probability density Moreover, this function satisfies a condition 



lim Uni (r) = 0. 



r-*-0 



(37) 



Next, the equation is transformed into dimensionless form. Introducing dimension-free en- 
ergy and potential, one obtains an equation, where each member has dimension that equals 

1/ {length f : 

Uni (r) = 0, (38) 



+ ^ ^mc^e {cni - V (r)) 



where e is combination of Hamiltonian parameters with energy dimension. For deep well 
and Coulomb potentials this parameter equals mc^. For Harmonic Oscillator potential 



V (r) = muj r /2 



(39) 



parameter e — hw is the best choice. It follows that parameter with length dimension for 
our problem is 



,2 ^ (^c) 
mc^e 



Taking new radial variable 



the equation converts to form 



d2 /(/ + !) ^ ^ (e„,-v(60) 



Unl (0 = 0. 



de e 

The discrete spectrum eigenf unctions of this equation form the orthonormal set: 

I K'l (0 «n/ (0 di = 8n',n- 

Jo 

The kinetic energy expectation value is 

1 r°° 

(to)„, = -^e (0 
The potential energy expectation value equals 



(40) 
(41) 

(42) 



(43) 



d^ _/(/ + !) 

de~^^ 



Unl (C) dC 



POO 

{^)nl = / <1 (0 V m Unl (0 di. 

Jo 



(44) 



(45) 
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Let us now compare bound eigenvalues spectra of two Hamiltonians - one with pseudo- 
relativistic kinetic energy operator, i.e. 



i-O — Lo 



to + V 



(46) 



and second - with relativistic kinetic energy operator 



h = t + V, 



(47) 



defined in f|T9l) . Let us mark the eigenvalue of the first Hamiltonian as e„;. The eigenvalue 
of the second Hamiltonian in first perturbation theory approximation equals 



{^)nl = ^nl + (Ah) 



nl ' 



(48) 



where, according to given above expressions, the eigenvalue correction due to kinetic energy 
relativistic expression is 



(Ah),, = (t - to),, = mc' 



1 + 2 (to),, /mc2 - 1 - (to),, /mc' 



(49) 



The first basic potential, which is going to be analyzed, is spherically symmetric well 
of infinite depth. This potential describes free particle, moving inside of sphere with no 
probability to find particle penetrated to outer region and equals: 



V r 



0, if r <d 
oo, if r > d 



(50) 



where d is sphere radius. The solution of corresponding Schrodinger equation (H2|) at v (bC,) = 
is 

«nK0=O/(0, (51) 

where ji (^) is the spherical Bessel function (j4|, § 10.1). The quantum number n = 1,2, ... 
is defined according to the boundary condition 



Uni (d/b) = 0. 



(52) 



This equation has an infinite number of solutions, defined by different values of 6 = 6„. For 
/ = 

Uno (0 = Nnosin^, (53) 
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nir. 



hence 

bn = d/ nil, Nno = a/27 
Thus, the eigenvalues of Schrodinger equation ( H2l) equal to: 

eno/mc^ = z {nnf /2, 



where z is the dimensionless well parameter 



he 



(54) 



(55) 



(56) 



d mc^ 

Potential equals zero in an inside region, hence the kinetic energy expectation value for this 
potential equals eigenvalue e„o- The perturbation due to relativistic corrections is 



1 + z (nvr)^ — 1 — z {uTc)'^ /2. 



This equals one percent of eigenvalue at 



zimrY ~ 0.083. 



(57) 



(58) 



For nucleon (mc^ ^ 1 GeV), moving in well, which radius equals 1 fm, the parameter 

^)^--^ (-) 

so this kinetic energy correction give a remarkable effect even for ground state [n = 1) 
of corresponding strong process model. After kinetic energy correction, the spectrum of 
spherical well will change in comparison with the well-known levels distribution, whose 
energies are proportional to n^. The distances between pairs of neighboring levels will be 
smaller than in spectrum without relativistic kinetic energy correction. 

Three-dimensional Harmonic Oscillator potential is next on the list. As previously de- 
fined, the parameters of the equation equal: 



e = hu, eb^ = {hcf /mc^, b"^ = h/muj. 



So, 



V (60 = eeV2, 



and the dimensions-free equation is 



rf2 / (/ + 1) 



de 



Unl (0 = 0. 



(60) 



(61) 



(62) 
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The eigenvalues of this equation are 



eni = huj{2n + I -1/2). (63) 



The eigenf unctions, normahzed by condition (H3|) are defined as in (j4!], § 22.6.18): 

uni (0 = N^i^^'e-^'^^LltT ie) . (64) 

Here normahzation constant 



^ 2"+'+i(n-l)! 
= v/^[2(n + /)-l]!!' ^^^^ 



Laguerre polynomial is defined as 



fc 



j=0 

where 



/3 \ _ r (/3 + 1) 
J 



(67) 



j!r(/3-j + i) 

is a binomial coefficient and F (/3) — Gamma function [T (/3 + 1) = /3r (/3) , F (1/2) = 
The potential expectation value equals the kinetic energy expectation value and both are 
defined as 

{^)ni /^c' = {^o)ni = eni/2mc^ = z{2n + l- 1/2) /2, (68) 
where dimensionless parameter 

ic \^ 
b mc^ J 

The Hamiltonian eigenvalue relativistic correction, defined in ( H9ll is: 



hio/mc^. (69) 



(Ah)„j /mc^ = VI + z (2n + / - 1/2) - 1 - z {2n + I - 1/2) /2. (70) 

As in square well case, this correction again will produce one percent of bound state energy 
at 

2(2n + /- 1/2) ^ 0.083. (71) 

As in earlier evaluation, at 2; ~ 0.04, this correction will be remarkable even in ground state 
of strong interaction models with such Harmonic Oscillator potential, i.e. at n = 1, / = 0. 
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The spectrum of corrected Hamiltonian will not consist of equidistant levels. The levels 
will be distributed closer than in the well-known Harmonic Oscillator spectrum with usual 
pseudo-relativistic kinetic energy operator. 

The last basic Hamiltonian is for Hydrogen-like atom with infinitely heavy nucleus. The 
Coulomb potential equals 



V (r) = —aZhc-, 



(72) 



where a ~ 1/137.036 is fine structure constant, Ze - nucleus charge and he ~ 0.197 keV 
nm. The length parameter of this problem equals Bohr radius: 

he 



aZmc^ ' 



the dimensionless problem parameter 



he 



a ■ me^ 



{aZf . 



After the new variable definition 



the dimensions-free equation looks as 



i = r/a, 



f_ _ + ^ 2 ^ 2eni 



Unl (0 = 0. 



The eigenvalue equals 



The corresponding eigenf unction is 



|], § 22.6.17): 



u 



nl (0 = iV.,p'+^e-''/^L^Y-i (P) 



(73) 



(74) 



(75) 



(76) 



(77) 



(78) 



where 



and normalization factor 



p = 2e/n 



in -I -I) 



(79) 



(80) 



(n + /)! 

The potential energy expectation value for this problem is two times larger than the eigen 
value 



n 



(81) 
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so the kinetic energy mean value is 




(82) 



The expression for eigenvalue correction now is: 



(Ah)„, /mc^ 



1 - z/2n^. 



(83) 



This equals one percent of energy eigenvalue when 



Z/n ^ 39.48. 



(84) 



It happens only at extremely hard conditions. It reaches maximum value for ground state 
only at large nucleus charge. This corresponds to the known character of relativistic effects, 
significant mainly for heavy atoms ([sl], § 11.7). 

IV. CONCLUSIONS 

Our consideration of relativistic corrections for basic quantum Hamiltonian has shown 
that the well-known momentum operator is not completely non-relativistic, as it has been 
widely believed. Pseudo-relativistic kinetic energy operator, present in Schrodinger equa- 
tion produces larger kinetic energy expectation values than the relativistic operator. It 
appears that relativistic kinetic energy correction can produce remarkable Schrodinger equa- 
tion spectrum transformations. The further consideration of Dirac equations transformation 
into Schrodinger form has shown that kinetic energy modification, known as weakly rela- 
tivistic potential correction, can not be useful for actual problems, because modern strong 
interaction realistic potentials have already all the possible potential members, allowed by 
symmetry consideration. 
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